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REPORT  SUMMARY 


The  purpose  of  this  research  project  was  to  evaluate 
an  optical  technique  for  the  determination  of  ion  densities  in 
a  high  altitude  chemical  release  of  barium.  The  method  involved 
the  detection  of  the  isotopically  shifted  hyper fine  line3  of 
barium  ions  whose  emission  was  less  likely  to  be  optically  thick 
due  to  their  minor  abundance  and  lower  relative  transition  prob¬ 
ability.  Thus,  emission  intensities  could  be  readily  related 
to  the  ion  population.  Simultaneous  detection  of  a  neutral  and 
an  ion  emission  line  could  also  be  accomplished  with  this  instru¬ 
ment  to  provide  data  regarding  the  ionization  conversion  process 
at  early  times. 

The  basic  instrument  is  a  Fabry-Perot  interferometer 
which  is  magnetos trie tively  scanned  offering  high  wavelength 
resolution  and  high  background  rejection.  The  instrument  and 
its  associated  equipment  were  government  furnished  through  the 
Defense  Atomic  Support  Agency.  Added  equipment  consisted  of 
narrow  band  optical  filters  for  the  bari- m  lines  of  interest. 

Data  were  obtained  for  the  releases  Apple  and  Dogwood 
but  sufficient  data  from  which  a  time  history  could  be  developed 
were  obtained  only  on  the  latter  event. 

Hyper fine  emission  features  were  observed  for  the  barium 
ion  line  at  4934  A  from  which  upper  excitation  level  populations 


were  dervied.  Calculations  were  made  to  determine  the  ratio 
between  total  ion  population  and  that  of  the  upper  energy  level 
giving  rise  to  the  4934  A  line.  Measured  intensities  are  then 
presented  as  barium  ion  column  densities  and  as  volume  densities 
by  consideration  of  the  pnotographically  derived  cloud  dimensions 
the  latter  results  agree  with  the  radar  measurements  of  electron 
density.  A  measure  of  the  neutral/ion  ratio  was  also  determined 
for  two  points  at  early  times.  The  results  thus  indicate  that 
the  proposed  technique  is  useful  in  the  determination  of  ion 
density  in  the  cloud  with  spatial  resolution  of  0.3  degrees. 


INTERFEROMETRIC  MEASUREMENTS  ON  SECEDE  I  BARIUM  RELEASES 


Robert  W.  Deuel  and  Robert  D.  Sears 
I IT  Research  Institute 
Chicago,  Illinois  60616 


I.  OBJECTIVES  OF  THE  EXPERIMENT 


A  high  resolution,  magnetos  trie tively  scanned 
Fabry  Perot  Interferometer  was  fielded  on  the  SECEDE  I  experi¬ 
ment  by  IIT  Research  Institute  personnel.  The  overall 
objective  of  the  IITRI  experiment  was  to  evaluate  an  optical 
measurement  technique  by  which  Barium  ion  densities  in  the 
cloud  may  be  determined.  The  high  spectral  resolution  (0.02A) 
of  the  instrument  allows  detection  of  isotope  shifted 
hyperfine  lines  of  Bal  and  Ball.  Because  of  the  lower  density 
of  the  minor  barium  isotopes  (135  and  137)  and  the  consequent 
longer  optical  path  in  the  cloud  for  the  hyperfine  lines, 
these  spectral  components  may  provide  an  "optically  thin" 
measure  of  Ba  neutral  or  ion  density  when  the  main  isotope 
(138)  lines  are  "optically  thick".  In  addition,  the  scanning 
capability  of  the  instrument  allows  simultaneous  intensity 
measurement  of  two  closely  spaced  lines  6498.78  and  6496. 9A 
which  are  scattered  by  Bal  and  Ball  respectively.  Thus  the 
Bal  to  Ball  conversion  process  may  be  followed  in  early  times. 
Finally,  the  high  spectral  resolution  of  the  instrument  allows 
detection  of  the  resonance  lines  even  against  a  high  background 
of  scattered  sunlight.  It  was  the  goal  of  this  part  of  the 
experiment  to  follow  the  Barium  ion  cloud  into  full  daylight. 


II.  DESCRIPTION  OF  THE  EXPERIMENT 

A.  Equipment 

The  interferometer  optical  unit  consists  of  a  single 
magnetos trictively  scanned  Fabry-Perot  etalon  preceded  by  a 
narrow  band  dielectric  filter  (on  the  order  of  0.1%  of  the  wave¬ 
length)  .  The  spacing  of.  the  interferometer  plates  is  about  5  mm 

0  0  o 

defining  a  free  spectral  range  of  0.24  A  at  4934  A  and  0.43  A 

O 

at  6498  A.  The  light  admitted  by  the  force  filter  thus 
contributes  about  20  orders  at  the  lower  wavelength  and  10  at 
the  higher  wavelength.  Selection  of  the  filters  and  features 
observed  prevented  unwanted  contamination  from  other  lines  in 
other  orders.  The  disadvantage  of  this  mismatch  occurs  at  late 
times  after  the  release  when  the  increasing  scattered  sunlight 
is  admitted  in  these  multiple  orders  thereby  decreasing  the 
optimum  signal  to  noise  ratio. 

The  optical  geometry  of  the  interferometer  unit 
defines  a  field  of  view  of  about  a  third  of  a  degree  with  a 
one  inch  collector  aperture.  Pointing  of  the  input  axis  is 
accomplished  by  a  two-axes  rotating  prism  mount  in  front  of 
the  collector  aperture.  Sighting  is  performed  through  a 
parallel  prism  and  telescope  arrangement.  The  fore  filter  for 
the  instrument  is  mounted  in  a  sliding  tray  which  has  provision 
for  two  separate  filters  allowing  sequential  observations  of 
two  features. 

The  scanning  interferometer  system  was  designed  to 
integrate  successive  wavelength  scans  in  a  coherent  manner  by 
incorporating  a  mercury  trigger  line  on  an  off  axis  path.  The 
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light  enters  through  a  fiber  optic  pipe  at  the  field-defining 
input  aperture  and  is  detected  through  another  fiber  optic  pipe 
at  the  corresponding  point  at  the  output  aperture.  Detection  of 
the  trigger  signal  thus  provides  an  accurate  wavelength  calibration 
for  subsequent  data  in  each  scan  of  the  etalon.  In  order  to 
prevent  contamination  of  the  main  detector  output  from  scattered 
mercury  light,  the  trigger  light  source  is  blocked  by  a  shutter 
after  the  line  i?j  detected. 

The  main  photomultiplier  output  observes  the  central 
fringe  at  the  output  aperture  which  corresponds  to  different 
wavelengths  as  the  interferometer  plates  are  magnetostrictively 
scanned.  To  enhance  detection  of  weak  signals  against  a  strong 
background,  the  output  signal  is  digitally  recorded  in  a  multi¬ 
channel  analyzer  where  each  pulse  represents  a  photon  interaction 
in  the  PM.  The  time  dependent  sweep  of  the  analyzer  and  the 
current  through  the  interferometer  bars  produces  a  direct 
correlation  between  channel  number  and  wavelength.  The  mercury 
trigger  signal  is  then  used  to  start  the  analyzer  sweep  at  the 
same  wavelength  point  for  coherent  addition  of  successive  scans. 
This  real  time  integration  and  display  of  the  channel  count 
permitted  integration  until  the  desired  spectral  feature  was 
observed  with  subsequent  readout  or  storage  of  the  data  to  provide 
maximum  time  resolution.  This  technique  was  used  during  the.  field 
tests  to  simultaneously  monitor  total  system  operation  and  to 
select  the  optimum  data  accumulation/readout  cycle. 

The  system  was  geneially  operational  prior  to  this  field 
III  RESEARCH  INSTITUTE 
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experiment  and  it  was  only  necessary  during  the  first  phase  of 
the  work  to  set  up  the  system  in  the  laboratory  for  check  out. 

New  coatings  were  required  for  the  interferometer  plates  to 
provide  the  required  spectral  resolution  throughout  the  range  of 
interest;  the  resolution  of  the  instrument  is  directly  related 
to  the  reflectivity  of  the  coating  at  the  desired  wavelength. 

By  this  stage  of  the  preparations,  it  was  decided  that  the 

O 

primary  features  to  be  observed  were  the  ion  lines  at  4934  A 

O  o 

and  6496  A  and  the  neutral  line  at  6498  A;  the  latter  two  lines 
would  be  observed  simultaneously  in  different  orders  through  the 
same  fore  filter.  It  was  thus  necessary  to  specify  high  (90%) 

O  O 

reflectivity  from  about  5000  A  to  6500  A  for  the  plates  of  the 
etalon.  The  requisite  narrow  band  fore  filters  were  also  ordered 
at  this  time. 

Successful  check  out  of  the  system  was  completed  prior 
to  shipment  of  the  equipment  to  Ramey  Air  Force  Base,  Puerto  Rico. 
The  site  chosen  to  provide  maximum  observational  capability  was  the 
Foreign  Broadcast  Information  Service  facility  near  Corozo  in 
southwestern  Puerto  Rico  near  one  of  the  photo triangulation 
sites.  The  equipment  arrived,  was  set  up  and  completely 
operational  during  the  week  before  the  planned  first  shot  of 
May  2. 


Ill  RESEARCH  INSTITUTE 
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B.  FIELD  OPERATION 


The  following  log  summarizes  activity  during  the  test 

series . 


May  2  Thur . 


May  4  Sat. 

May  5  Sun . 

May  7  Tues . 

May  8  Wed. 

May  12  Sun . 


Operational  02:30;  attempted  calibration 
aborted  due  to  condensation  on  interior 

O  O 

optical  components;  4934  A  and  6498  A 
filters  installed.  Launch  0512. 

O 

One  data  record  obtained  for  4934  A, 
integrated  from  T+20  to  T+240  seconds. 
Operational  02:4.5;  mission  scrubbed  at 
T-10  minutes,  system  checks  with  sodium 
lamp. 

Operational  02:45;  mission  scrubbed  03:45. 
Operational  02:30;  launch  05:00  (-12.5°). 
Operate  until  05:55.  Analysis  indicates 
no  data  obtained  this  shot. 

Operational  03:00;  mission  scrubbed  03:45; 
azimuth  and  elevation  corrections  for 
Corozo  noted  with  respect  to  Parguera. 
Operational  03:30;  launch  05:06;  second 
stage  observed,  no  release,  and  no  data. 
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May  13  Mon.  Operational  03:20;  launch  05:08;  operate 

until  06:14;  8  data  records  for  4934  A 

O 

to  T+24  minutes  and  for  6498  A,  three 
records  from  T+O  to  T+150  seconds. 
Calibration  for  resolution  performed 
after  completion  of  test  data. 


III. 


DESCRIPTION  OF  THE  DATA 


A.  Acquisition 

The  release  clouds  were  successfully  acquired  on  the 
telescopic  viewfinder  of  the  interferometer  on  Haven  Barter  I 
and  IV.  On  Haven  Barter  II,  the  release  cloud  faded  from  view 
before  the  interferometer  field  of  view  could  be  adjusted  pro¬ 
perly  in  azimuth  and  elevation.  No  release  took  place  on 
Haven  Barter  III. 

During  the  first  event,  a  35  mm  camera  was  attached  to 
the  tracking  telescope  to  provide  specific  data  regarding 
pointing  direction  for  correlation  to  the  line  emission  measure¬ 
ments.  Visual  sighting  for  manual  tracking  was  then  accomplished 
through  the  reflex  optics  of  the  camera.  It  was  immediately 
apparent  during  the  first  event  that  the  extent  of  the  barium 
cloud  and  the  relatively  low  contrast  obtained  at  the  camera 
view  finder  made  tracking  extremely  difficult.  Thus,  the 
paucity  of  data  during  the  first  event  is  at  least  partly 
ascribed  to  poor  tracking  of  the  cloud.  The  camera  was  removed 
from  the  sighting  telescope  for  the  subsequent  events  where  the 
pointing  direction  was  indicated  on  sketches  of  the  cloud 
structure.  Tracking  problems  remained,  however,  because  the 
telescope's  field  of  view  was  smaller  than  the  cloud  size 
except  at  very  early  times  and  also  because  of  the  inability 
to  visually  locate  the  cloud  at  late  times.  These  problems 
very  likely  contributed  to  the  variance  of  data  during  the 
observation  as  well  as  the  termination  of  data  when  the  back- 
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ground  increased. 

Observations  were  made  during  the  second  event  but 
analysis  of  the  data  revealed  no  evidence  of  line  emissions. 
Although  operational  for  the  third  event,  the  apparent  mal¬ 
function  of  the  payload  precluded  data  acquisition. 

The  fourth  event  produced  eight  sets  of  intensity 

O 

measurements  for  the  4934  A  ion  line  from  T+3  minutes  to  T+23 

O 

minutes  and  two  sets  of  intensity  measurements  for  the  6498  A 

O 

neutral  and  6496  A  barium  ion  line  from  T+30  seconds  to  T+180 
seconds . 


B.  Hyper fine  Structure  and  Line  Shape 

The  interferometer  scanning  interval  was  set  to  cover 
slightly  more  than  two  orders  for  both  the  4934  A  line  and  the 
6496.9  and  6498.78  A  lines.  Hence,  placement  of  the  lines  with 
respect  to  the  free  spectral  range  resulted  in  reproduction  of 
three  main  peaks  for  the  4934  A  line,  and  two  sets  of  peaks  of 
the  6496.9  and  6498.78  A  lines.  The  spectral  intensity  which 
is  expressed  as  a  count  rate  versus  wavelength  in  terms  of  the 
counting  channel  assignment  is  illustrated  in  Figure  1.  Because 
the  wavelength  scan  is  somewhat  non-linear  in  relation  to  the 
driving  current  and  hence  non-linear  in  channel  assignment  (time), 
the  third  order  peak  at  the  right  side  of  the  4934  A  display 
is  wider  than  the  first  and  second  orders.  The  same  effect  may 
be  noted  on  the  6496.8  A  plot. 

Instrumental  characteristics  determine  the  ability  to 
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N  (counts /sec. 
(Rayleighs)  -  N  x 


resolve  the  hyperfine  structure  on  the  4934  A  line.  The 
hyperfine  lines  are  located  at  separations  of  ±  0.026A, 

+  0.034  A,  and  +  0.038A  from  the  main  emission  peak.  The 
instrumental  resolution  for  the  4934  A  and  for  the  6498  A 
region  is  determined  by  the  free  spectral  range  and  by  the 
finesse  of  the  etalon  plates.  These  quantities  are  listed 
in  the  following  table  together  with  the  observed  line 
characteristics . 

TABLE  1 

WAVELENGTH  FREE  SPECTRAL  RANGE  FINESSE  RESOLUTION  OBSERVED 

LINE  WIDTH 

4934  A  0.24  A  17  0.014  A  0.04  A 

6498.78  A  0.43  A  13  0.033  A  0.05  A 

The  observed  line  widths  for  the  main  emission  peaks  were  much 
larger  than  expected  insofar  as  the  doppler  width  of  0.008A 
(corresponding  to  an  ionospheric  temperature  of  800°K)  would 
place  the  natural  line  within  one  or  two  time-spectral  resolution 
channels.  The  instrumental  characteristics  were  carefully 
examined  in  order  to  ascertain  positively  that  the  line  width 
was  not  primarily  an  internal  effect.  Lack  of  any  clear-cut 
instrumentation  line  broadening  effects  leads  to  the  conclusion 
that  the  observed  line  breadth  is  real  and  that  it  is  an  effect 
of  the  optical  depth  of  the  emitting  medium.  Comparison  of  the 
line  widths  with  an  idealized  Gaussian  profile  for  a  800°K 
doppler  broadened  line  indicates  that  the  true  peak  emission 
intensity  is  a  factor  of  at  least  2  to  3  higher  than  the  observed 


10 


peak.  As  will  be  seen,  this  conclusion  conforms  with  that  made 
by  comparison  of  the  main  line  observed  intensity  values  with 
the  intensity  of  the  hyperfine  structure  lines. 

The  instrumental  resolution  is  adequate  to  observe  the 
0.026  A  and  the  0.034/0.038  A  lines.  As  seen  in  Figure  1,  these 
lines  are  of  much  lower  intensity  than  the  main  line,  however, 
they  appear  consistently  in  all  of  the  4934  A  data  in  the  same 
wavelength  channels.  It  is  possible  to  determine  the  intensity 
of  the  hyperfine  emissions  by  subtraction  of  the  main  line 
profile  and  background  count  rates  from  the  peak  count  rate  in 
the  hyperfine  line  channels.  Because  the  hyperfine  lines  are 
expected  to  be  optically  thin,  and  they  lie  in  at  most  two 
wavelength  channels,  the  single  channel  peak  value  is  a  reasonable 
estimate  of  the  total  hyperfine  line  intensity. 

The  results  for  the  fourth  event  (6804  Dogwood) ,  shown 

O 

on  Figure  2,  indicate  a  roughly  constant  level  for  the  4934  A 
primary  peak  (Ba  138)  throughout  the  observation  period  from 
T+3  to  T+23  minutes.  It  is  noted  that  the  largest  deviation 
from  the  constant  level  at  T+5.5  minutes  corresponds  to  an 
observation  near  the  edge  of  the  cloud  while  all  other  measure¬ 
ments  correspond  to  positions  near  the  center  of  the  ion  cloud. 

For  comparison,  the  isotope  line  intensity  is  shown  normalized 
to  account  for  the  difference  in  relative  abundance  of  the 
respective  species  (135/137  to  138)  and  for  the  difference  in 
transition  probability.  The  intensity  of  the  barium  ion  isotope 
line  appears  to  decrease  during  the  first  twelve  minutes  with 
the  subsequent  two  data  points  at  T+20  and  T-t-22  minutes  indicating 
an  increase.  Although  the  error  limits  are  on  the  order  of  a 


factor  of  two  or  three,  this  interpretation  of  the  time  history 
is  consistent  with  the  data  of  the  primary  138  isotope  which 
indicates  optical  thickness  at  times  less  than  T+10  minutes. 


Figure  1  shows  average  data  for  both  the  6498/6  A  and 
4934  A  lines;  the  upper  curve  is  the  integration  of  the  three 
sets  of  data  totaling  100  scans  (at  one  second/scan)  for  the 
6498/6  A  neutral  and  ion  lines;  the  lower  curve  is  the  integrated 
data  for  five  of  the  eight  sets  of  measurements  for  the  4934  A 
ion  line  totaling  300  scans.  The  primary  feature  observed  in 
the  first  case  is  the  neutral  barium  emission  at  6498.76  A  seen 
in  two  successive  orders.  At  approximately  ten  and  twenty 
channels  to  the  left  of  this  feature  for  the  first  and  second 
orders  respectively,  there  is  observed  a  minor  increase  in 
intensity  in  the  appropriate  channel  for  the  ion  line  at  6496.90  A. 
Assuming  this  identification,  the  average  background  is  subtracted 
from  these  channels  to  provide  ion  intensity  data  corresponding 
to  the  readily  determined  neutral  line  intensity.  In  Figure  2, 
two  points  are  plotted  for  each  of  these  lines  showing  the  decay 
of  the  neutral  and  the  growth  of  the  ion  emission.  A  third  point 
for  each  line  at  T+2  minutes  is  in  the  noise  level  negating 
comparison.  However,  the  first  two  points  may  be  used  to  obtain 
decay  and  growth  rates  for  the  neutral  and  ion  emission  respectively. 
An  exponential  fit  to  these  points  provides  a  1/e  decay  time 
for  the  6498  A  neutral  barium  line  of  about  24  seconds  and  an 
e-folding  increase  for  the  6496  A  barium  ion  line  of  about 
60  +  30  seconds.  It  is  noted  here  that  the  figure  shows  count 
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rates  per  second  which  are  related  to  intensity  by  the  trans¬ 
mission  of  the  fore  filter  at  the  wavelengths  of  the  spectral 
lines.  The  varying  transmission  results  in  an  intensity  ratio 
of  two-thirds  the  count  rate  ratio.  Thus, 

TABLE  2 

Time  6498  6496 

Intensity  (Neutral/Ion) 
Ratio 

T+42  sec .  12 

T+72  sec.  4 

When  the  signal  at  6498  A  was  lost  in  the  noise  level 
at  T+2  minutes,  the  interferometer  was  switched  to  the  4934 
filter.  Subsequent  observations  were  made  at  6498  A  but  no 
further  data  were  obtained. 
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IV. 


DETERMINATION  OF  ION  POPULATION 


The  direct  data  obtained  in  this  experiment  in  counts 
per  second  per  unit  wavelength  correspond  to  the  input  intensity 
in  photons  per  second  within  the  same  wavelength  interval.  The 
photon  counting  technique  is  better  than  90%  efficient  for  the 
count  rates  encountered  here.  Considering  overall  optical 
efficiency  of  the  instrument,  the  quantum  efficiency  of  the 
photomultiplier,  the  field  of  view,  and  the  collector  aperture, 
the  following  relation  is  obtained. 


^source 


shotons 


sec  cm 


-  1.82  x  108  x  N 


It  is  noted  here  that  the  derivation  of  this  relation  involved  the 
correlation  of  the  count  rate  at  late  times  (T+30  to  40  minutes 
when  the  signal  is  due  to  scattered  sunlight)  with  published  data 
for  scattered  sunlight  versus  solar  angle.  The  empirical  pro¬ 
portionality  constant  thus  contains  a  lumped  parameter  including 
lens  and  filter  transmissions,  mirror  and  prism  reflectivities, 
interferometer  transmission,  and  PM  quantum  efficiency.  This 
pseudo-calibration  procedure  likely  provides  the  most  realistic 
conversion  due  to  the  inherent  inclusion  of  the  total  system 
function  but  it  is  recognized  that  greater  precision  with  proper 
calibration  is  highly  desirable.  This  point  should  be  remembered 
in  the  interpretation  of  the  derived  data. 


A.  Intensity  versus  Number  Density  (Excited  State] 


A  linear  relation  between  observed  intensity  and  the 

number  density  of  the  upper  level  giving  rise  to  the  observed 
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radiation  holds  for  the  region  where  the  line  is  said  to  be 
optically  thin.  The  addition  of  more  emitters  along  the  line- 
of-sight  path  does  not  then  produce  a  corresponding  increase  in 
intensity  due  to  the  occurrence  of  self-absorption.  The  relation 
between  observed  intensity  and  line  density  (volume  density  x 
path  length)  can  be  calculated  for  a  pure  Doppler  broadened  line 

O 

(Plass  1958)  and  the  result  is  shown  in  Figure  3  for  the  Ball  4934  A 

line.  The  linear  relation  holds  for  (NL)  less  than  10^  while  beyond 

13  1/2 

10  the  intensity  is  prof  rtional  to  (In  NL)  ;  the  latter  region 

is  where  the  line  is  said  to  be  optically  thick. 

Observation  of  the  hyperfine  lines  of  the  Ball  4934  A  line 

should  provide  a  source  whose  intensity  is  about  a  factor  of  20  less 

than  the  primary  line  andwhich  is  thus  more  likely  to  be  optically 

thin.  From  the  electron  density  data  at  times  corresponding  to  the 

6  7 

optical  data,  the  ion  number  density  is  between  10  and  10  with  an 

ion  cloud  radius  between  2  and  3  km  (Rosenberg) ;  the  resulting  ion 

11  12 

line  density  is  5  x  10  to  5  x  10  .  Considering  the  factor  of 

twenty  reduction  in  relative  emission  for  the  isotope  line,  its 
intensity  should  fall  in  the  region  shown  on  the  lower  curve  in 
the  figure.  The  derived  values  for  ion  line  density  may  later 
be  examined  on  this  figure  to  determine  if  corrections  are 
necessary  for  optical  thickness. 

The  proportion  of  excited  state  ions  for  the  Ball  4934 
isotope  transition  to  the  total  barium  ions  in  this  state  is 
1/21.9.  The  line  density  is  then  found  from  the  isotope  emission 
intensity  as  follows: 
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Figure  3 

O 

Intensity  versus  Line  Density  for  Ba  II  4934  A.  Emiss'.on 


2 

where  is  the  number  density  of  barium  ions  in  the  P^/2  ^-eve^> 
L  is  the  optical  path  length,  and  A41  is  the  Einstein  transition 

O 

probability  for  4934  A  emission.  In  order  to  determine  the 
total  ion  population, it  is  necessary  to  determine  the  relative 
populations  of  barium  ions  in  all  possible  states  to  find  the 
ratio  of  N^/N^tai*  The  following  section  will  discuss  the 
excitation  processes  and  provide  calculations  of  the  ratio  of 
total  ions  to  ions  in  the  upper  level  of  the  4934  A  line. 

B.  Population  of  Ball  States 

Both  radiative  and  coilisional  excitation  processes 
will  be  considered  to  determine  the  relative  excitation  rates. 


B.l  Coilisional  Excitation 


Considering  two  Maxwell -Boltzmann  distributions,  one 
of  Ball  (ground  state  number  density  »  N^)  and  one  for  an  electron 
and  proceeding  in  a  manner  similar  to  that  indicated  in  Chapman 
and  Cowling  (I960) ,  one  obtains  the  following  expression  for 
the  rate  of  coilisional  excitation  of  Ball: 


dN 


14 

3t“ 


2 

N,N 
fr  1  e 


2  kT 
— 3F" 
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1/2 

-)  (T 


*  2 
m  g  . 


ae 


(1  +  — 2'  kT")  exp 


*  2 
-®  «o 

TW 


*  ml  me 

whore  m  =  —  5~  >  go 

I  e 

ponding  to  the  minimum  energy 


is  the  relative  velocity  corres- 
bolow  which  no  excitation  occurs, 
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and  (T  is  the  atom-electron  collision  cross  section.  The 
ae 

derivation  of  this  expression  assumes  that  the  same  temperature 

can  be  applied  to  each  distribution;  the  case  for  unequal 

temperatures  is  more  complex  and  will  not  be  considered  here. 

•k 

With  ‘.ne  appropriate  cross  section  (  CTl^)  and  energy  ,  m  ^o  \ 

ae  v  2  / 

equation  (3)  gives  the  rate  of  excitation  to  a  given  level. 
Consider  the  excitation  process: 


Ba+  (2S,  /0)  +  e  -»-Ba+(2P1  /9)  E  =  20261.5  cm" 


The  barium  ion  is  essentially  stationary  compared  with  the  electron, 

so  that  gQ  must  be  a  limiting  electron  velocity.  In  the  center 

of  mass  system  used  above,  the  electron  can  give  up  only  one-half 

k 

of  its  energy  maximum.  Hence  m  &o  must  equal  twice  the 

- 

excitation  energy  or; 


*  2 


m  go  _  2  x  20261.5  x  6.6256  x  3.00  x  10~17  ,= 


2  kT 


1.3305  x  800  x  10 


=  72.9 


6  -3 

The  electron  density,  Ng,  is  **10  cm  .  Thus  equation  (3) 
becomes 


dN14 


=  Nx  (Tae  x  1  106  (7.22  x  1014)  1/2  x  72.9  x  10' 


=  N,  <T  x  1.2  x  10”16  sec.”1  cm”2 

i.  dc 
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•1  f)  ? 

Since  Q~afi  is  on  the  order  of  10  cm  ,  equation  (5a)  reduces 
to: 


<”14 

dt 


-VI 

Nx  x  10 


sec , 


-1 


B.2  Radiative  Excitation 


(5b) 


2 

Population  of  the  P-^/2  state  from  the  ground  state 
can  also  occur  by  the  absorption  of  radiation  according  to  the 
relation: 


dN,, 

crr-=  Ni  Bu  *u  <»)  (6a) 

where  N-^  is  the  number  density  in  the  lower  state 

Bij  is  the  Einstein  transition  probability  for  absorption  and 
j>l^(U)  is  the  spectral  radiation  density  at  the  wavelength 
of  the  transition. 


The  transition  probability  for  absorption  (Bij)  is  proportional 
to  the  transition  probability  for  emission  (Aij): 


B 


12 


8  h  )} 


kry,  =  6.55  X  10 
i  1. 
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cm 


erS  sec ^ 


The  excitation  sour ;e  of  radiant  energy  is  sunlight  which  has  a 

o 

spectral  density  at  5000  A  of: 


0.20  x  10~4  x  107 

V* 


cm 


sec 
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_  90  _  •} 

=  5.4  x  10  ergs.  sec.  cm 

The  radiative  excitation  rate  is  then: 

Nx  x  6.55  x  1019  x  5.4  x  10"20  sec"1 

x  3.55  sec"1  (6b) 

These  rather  detailed  calculations  for  two  kinds  of  excitation 
processes  provide  the  basis  for  the  following  discussion  of 
the  order  of  magnetude  for  all  possible  kinds  of  excitation  and 
de-excitation  processes,  which  may  affect  the  relative  populations 
of  the  states  for  the  barium  ion. 

B.3  Other  Excitation/De-excitation  Processes 

All  of  the  events  that  can  populate  or  depopulate 
excited  states  must  now  be  considered  if  we  are  to  relate  an 
excited  state  ion  concentration  to  the  total  ion  concentration. 

We  will  let  Nj  be  the  concentration  of  ions  in  the  jth  excited 

O 

state  (in  our  case  the  upper  state  of  the  4934  A  transition). 

N-j^  is  the  ground  state  population,  the  population  of  the 
2 

metastable  D  levels,  the  combined  population  of  all  states 

above  Nj  ,  is  the  electron  density,  and  OC,  p,  and  are  the  various 

excitation  and  de-excitation  rates.  All  of  the  results  of  this 

discussion  are  presented  in  Table  3. 

For  the  case  of  collisional  excitation,  the  OCwas  evaluated 

-16  2 

in  the  preceding  section  using  a  cross-section  of  /-w10  cm  . 
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Table  3  shows  the  value  of  10‘32  n  obtained  (event  type  1). 

For  event  type  2,  ion  and/or  atom  collision  excitation  events 

can  be  considered  approximately  equal;  the  cross  section  will  not 

be  different  by  more  than  one  or  two  orders  of  magnitude  from 

that  for  atom-electron.  Therefore,  it  is  reasonable  to  assume 

that  event  (1)  event  (2).  Events  (3)  and  (4)  refer  to 

2 

collisional  excitation  from  the  D  states;  calculations  indicate 
that  even  if  N^>N^,  an  upper  limit  to  this  inequality  would 

g 

be  ~10  .  Calculations  for  (5)  and  (6) ,  based  on  reasonable 
values  for  solar  irradiance  and  the  Einstein  B  coefficients 
provide  estimates  of  unity  times  the  lower  state  number  density. 

It  is  easily  seen  then  that  radiative  excitation  completely 
dominates  and  collisional  processes  can  be  neglected  in  the 
excitation  scheme. 

The  remaining  event  types  refer  to  the  possible  loss 
2 

modes  of  ions  from  the  P^/2  leve^*  Event  (7)  represents  the 

loss  from  state  j  due  to  electron  collisional  excitation  to 

higher  states  N,  .  Since  it  is  unlikely  that  the  j  state  population 

could  ever  become  as  large  as  the  ground  state  population,  it  is 

apparent  that  (7)  < (2) .  Event  (8)  is  less  than  event  (7)  because 

of  the  lower  relative  value  of  the  cross  section  for  these 

2 

respective  events.  Event  (9)  may  be  about  10  larger  than 
event  (2),  because  of  the  possibility  of  resonance  radiation 
quenching,  but  certainly  no  more  than  this.  The  direct  radiation 
loss  (10)  depends  simply  upon  the  radiative  lifetime,  hero,  assumed 
to  be  10  .  Events  (11)  and  (12)  represent  the  radiation  absorption 

loss  and  recombination  loss;  both  are  of  the  order  N j .  Thus 
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TABLE  3 


ORDER  OF  MAGNITUDE  COMPARISON  OF  COLLISIONAL  AND  RADIATIVE  EVENTS 
CONTROLLING  THE  POPULATION  OF  AN  EXCITED  ENERGY  STATE  j 


Type  of  Event 


No.  of  Events/  Quantitative  Estimates 

2 

cm  -sec  and  Comparisons 


(1)  Electron  Collision  ^i^l^e 

Excitation  from  the  J 

Ground  State 

(2)  Heavy  Particle  °S^l^i 

Collision  Excitation  J 

from  the  Ground  State 

(3)  Electron  Collision 
Excitation  from  the  kt“ 
state 

(4)  Heavy  Particle  Collision  °^NkN4 
Excitation  from  the  kfcb  J 
state 

(5)  Radiative  Excitation  from  f 

Ground  State  J 

(6)  Radiative  Excitation  from  f-, 

kth  state  3 


N 


1 


~(D 

10“24  Nx 

~(3) 


Sc 


(7) 

Electron  Colli s ion  loss  to 
higher  states  1 

Vj*. 

(7)  <  (2) 

(8) 

Super  elastic  electron 
Collisions 

(8)  <  (7) 

(9) 

Deactivitation  Collisions 
with  Heavy  Particles 

(9)~102r2) 

(10) 

Direct  Radiation  Loss 

10*  Nj 

(U) 

Raalation  Absorption  Loss 
to  Higher  States  1 

s 

(12) 

Recombination 

%  NjNe 

NJ 

radiative  processes  dominate  de-excitation  as  well  as  excitation 
processes . 

B.4  Reaction  Scheme 


In  order  to  find  the  total  ion  population  from  the 

intensity  of  one  ion  transition,  it  is  necessary  to  obtain  the 

2  2  2 

relation  between  the  populations  of  the  P,  D  and  the 

ground  state,  assuming  most  of  the  bar5.um  ions  are  in  these 

2 

states.  The  excitation  rate  of  the  levels  higher  than  P  as 

2 

shown  in  Figure  4,  must  be  Jbout  l/20th  that  of  the  P^/2  state 
due  to  the  diminished  spectral  radiance.  Hence  these  higher 
levels  will  be  ignored. 

The  steady  state  population  equations  are  then; 


N4A41  +  N5A51  "  N1B14^14  +  N1B15A5  ^ 

N4A42  +  N5A52  “  N2B24-^ 24  +  N2B25  ^25  (8) 

n5A53  “  N3B35^35  ^ 

NlB14fl4  +  N2B24^24  "  N4A41  +  N4A42  <10> 

N1B15  ^  15  +  N2B25^ 25  "  N5A51  +  N5A52 


Stimulated  emission  can  be  neglected  here  because  the  rate  relative 

8  -1 

to  that  for  spantaneous  emission  is  3.55/10  sec  .  The  N's  re¬ 
present  the  number  densities  of  the  respective  levels  shown  in 
Figure  4  and  Aij  is  the  transition  probability  for  emission. 
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Figure  4  Lower  Energy  Levels  of  Ball 


Transitions  of  direct  interest 

N's  represent  the  number  densities  in  each  state 


Gallagher  (1967)  lists  the  following  transition 


probabilities: 


TABLE  4 


Wavelength 

Designation 

Absolute 

Transition 

Probability 

4554.0 

A51 

1.18  x  108 

5853.7 

A52 

.048 

6141.7 

A53 

.37 

4934.1 

A41 

.95 

6496.9 

A42 

.33 

The  values  are  determined  from  the  following 
equation  (6a)  and  the  jy ,  (solar  spectral  radiation 
density) ,  are  obtained  from  the  compilation  by 
H.  H.  Malitson,  NASA,  Goddard  Space  Flight  Center. 
The  simultaneous  solution  for  equations  (7)  through 
(11)  then  provides  the  required  population  ratios. 
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Nx  =  2.14  x  108  Na 
N2  =  0.122  x  108  N4 
N3  =  3.11  x  108  N4 
N5  »  1.34  N4 

or  NTotal  *  5'37  x  1q8  n4 


(12) 


NOTE:  A  calculation  was  made  to  determine  the  relative 

2 

efficiency  of  radiative  excitation  of  the  P^/2  *-evel 
from  solar  radiation  versus  optical  pumping  from 
radiation  in  the  wings  of  the  Ba  136/138  unshifted  line. 
The  result  indicates  that  solar  excitation  would  account 
for  47%  of  the  total.  A  correction  for  this  effect  is 
included  in  the  reduction  of  intensity  to  ion  density. 
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C .  Barium  Line  Density  and  Volume  Density 

Combining  equations  1,2,  and  12  provides  the  relation 
between  count  rate  and  ion  line  density. 

vr  iOnS  _  r\  n  o  /\10  cts 

N  — T-  •  2'23  *  10  *  N  sec  (13- 

cm  v 

t 

Referring  again  to  Figure  3  and  recalling  that  the  original  data 

relates  to  the  intensity  of  the  hyperfine  isotope  line,  the 

effective  ion  line  density  below  which  the  radiation  may  be  said 

19  -2 

to  be  optically  thick  is  about  2  x  10  cm  All  of  the  data 
obtained  for  the  Ball  4934  A  line  lie  below  this  value  and  thus 
conform  to  the  linear  relation  to  number  density  of  the  emitting 
specie.  The  results  for  the  barium  ion  line  density  as  a  function 
cf  time  for  the  Test  6804  (Dogwood)  are  presented  in  Figure  5. 
Based  upon  the  preliminary  data  of  Rosenberg,  conversion  may  be 
made  to  volume  density  by  using  the  cloud's  optical  Gaussian 
halfwidth.  These  data  are  shown  in  Figure  6  along  with  the  volume 
electron  density  measurements  by  Raytheon.  The  large  scatter  of 
the  data  points  may  be  ascribed  in  part  to  the  variability  in 
pointing  with  respect  to  the  cloud  center  as  discussed  in  a  pro¬ 
ceeding  section.  A  linear  least  squares  fit  to  these  points  is 
shown  on  the  figure  and  indicates  a  slope  similar  to  that  for 
the  electron  density  change. 
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Lectron  and  Ion  Time  Dependence 


V, 


SUMMARY  AND  CONCLUSIONS 


Measurements  of  the  intensity  of  resonant  scattered 
emission  features  from  Barium  releases  in  the  upper  atmosphere 
have  been  successfully  related  to  columnar  ion  density  in  the 
Barium  cloud.  Through  adoption  of  photographic  estimates  of  the 
Gaussian  half  width  of  the  released  cloud,  a  peak  volume  density 
of  Barium  ions  has  been  derived.  Use  of  very  high  resolution 
interferometric  instrumentation  has  allowed  measurement  of  the 
isotope  shifted  hyperfine  structure  of  at  least  one  of  the 
lines  (Ball,  4934  A).  In  addition,  by  combining  line  measurements 
in  different  orders,  simultaneous  determination  of  Bal  and  Ball 
density  in  the  early  stages  of  cloud  development  was  shown  to  be 
practicable. 

General  agreement  of  che  temporal  behavior  of  the  peak 
Barium  ion  volume  density  values  derived  herein  with  the  electron 
density  data  which  were  obtained  from  high  frequency  radar  has 
been  demonstrated.  Although  absolute  value  of  peak  volume 
density  of  Barium  ions  is  a  factor  of  about  two  higher  than  the 
peak  electron  density,  this  difference  is  within  the  error  limits 
imposed  by  two  sources:  error  in  photographic  determination  of 
cloud  size,  and  a  systematic  error  source  in  the  calibration  which 
could  not  be  corrected  in  the  field.  Based  upon  these  results,  it 
is  concluded  that  the  optical  measurements  of  Barium  Ion  density 
will  provide  quantitatively  accurate  and  useful  results  in 
subsequent  experiments. 
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the  ratio  between  total  ion  population  and  that  of  the  upper  energy 
level  giving  rise  to  the  4934  A  line.  Measured  intensities  are 
presented  as  barium  ion  column  densities  and  as  volume  densities  by 
consideration  of  the  photographically  derived  cloud  dimensions.  A 
measure  of  the  neutral/ion  ratio  was  also  determined  for  two  points 
at  earlv  times. 
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